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THE  MQTIOH  OP  A PL 04. T IS G SPHERE  IE  SURFACE  RATES 

by 

R.  C.  MaoCaay 


Introduction? 

In  a reoent  report the  author  has  oonsidered  sane  approximate 
solutions  to  the  problem  of  a freely  floating  oiroular  oylinder  in  surfaoe 
gravity  waves,  ^he  objeet  of  the  present  report  is  the  treatment  of  the  analo- 
gous problem  for  a semi-immersed  sphere. 

The  general  problem  of  floating  obstaoles  has  been  treated  in  detail 
by  Frtts  JohnC^).  The  results,  however,  while  establishing  the  uniqueness 
and  existenoe  of  solutions,  are  oast  in  terms  of  integral  equations  whioh  in 
the  three  dimensional  oase  posssess  infinite  kernels.  This  being  the  case,  at- 
tempts to  obtain  numerical  solutions  along  these  lines  7oulu  to  prs3cn. 

formidable  problems. 

An  explicit  volution  in  the  oase  of  a semi-iasaersed  sphere  is  obtained 
in  this  report  by  a generalisation  of  the  method  of  separation  of  variables. 

The  method  is  based  os  the  procedure  used  by  Ursell'®^  in  the  case  of  the 
oiroular  oylinder.  In  order  to  satisfy  all  boundary  conditions  and  obtain  the 
oorreot  behavior  at  infinity,  the  velocity  potential  is  taken  as  a superposi- 
tion of  that  due  to  a point  ? cures  and  an  infinite  series  of  non- orthogonal 
harmonic  funotions.  To  obtain  the  motion  of  a freely  floating  sphere  we  first 
solve  the  problem  of  determining  the  motion  prc'dv.o**!  by  forced  oscillation.  A 
oorreotion  term  is  then  added  to  aooount  for  the  incident  wave,  and  a solu- 
tion is  obtained  which  contains,  as  unknowns,  the  amplitude  and  phase  of  the 
motion  of  the  sphere.  *hese  quantities  may  be  fixed  by  means  of  tha  equations 
of  motion. 

Sate  remarks  on  the  convergence  of  the  series  are  made  and  it  is 
established  that  the  series  in  the  feroed  oscillation  problem  converges, 
provided  the  frequenoy  of  motion  is  sufficiently  small,  the  proof  yielding 
at  the  same  time  an  alternative  method  for  the  calculation  of  the  coefficients 
in  the  expansion. 

Some  numerical  work  is  underway,  but  due  to  she  limited  computing 
facilities  available,  is  progressing  slowly.  In  view  of  the  suooess  of  the 
methods  in  the  report  on  the  oyilnaer , iv  would  he  exp; c ted  that  the  tech- 
nique should  lead  to  fairly  aoourate  results. 

1.  Formulation  of  the  Problemt 


We  suppose  an  incompressible,  non-visoc-i  fluid,  of  infinite  depth, 
to  fill  the  region  y>o  in  its  undisturbed  state.  *he  x-z  plane  is  then 
to  oednoide  with  the  free  surfaoe.  We  suppose  further  that  in  the  absenoe  of 
any  obstaoles  the  motion  of  the  fluid  is  that  of  a two-dimensional,  sinusoidal, 

* lumbers  in  parentheses  ref«ir  to  References  at  end  of  paper. 
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gravity  m?t  of  frequency  <r~  . If  now  & rigid  sphers,  of  radius,  a,  is  intro- 
duood  fnto  the  free  curfaoe,  it  will  koto  under  the  totion  of  the  inoaning 
wares.  In  its  equilibrium  position  the  sphere  is  to  t*  urai-iassersed  with  its 
oenter  at  the  origin.  Ice  notions  are  assuned  to  be  snail  and  periodio  wivh 
frequency  o~  • 


A.  o replete  determination  of  the  notion  is  possible  if  one  suooeeds 
in  finding  tL&  velocity  potential,  4 (x,y,*,t)  , whioh  is  a function, 
harmonic  in  the  domain  of  the  fluid  and  satisfying  oertain  boundary  conditions. 
Xho  snrfaos  elevation  abore  y s 0,  r{  (x,x,t),  and  the  pressure  p (x,y,£,t) 
are  obtained  from  4 by  neans  of  the  equations. 


L d* 


- 


The  assumption  of  sinusoidal  tine  variation  neans  we  nay  take, 
4 (x,y,s,t)  - He  (W(x,y,s)  e-1*-*) 
where,  within  the  framework  of  small  notions  W satisfies. 


*77 


* 0 in  fluid 


Wy*f  * 0 on  free  surface,  I = <J*z/g 
In  addition,  if  tq  is  the  normal  velooity  of  the  sphere, 

¥ = on  -2  . v2  ^ .2  . .2 


(3.1) 

(1.2) 

(1.3) 
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To  deal  with  ‘the  problez  of  the  freely  floating  sphere,  we  have 
first  to  solve  an  auxiliary  problem  in  whioh  the  sphere  moves  with  a forced 
▼ertioal  osoillation  of  froquenoy,  <r  , in  a fluid  initially  at  rest.  Intro- 
ducing spherioal  oo-ordinates  ), 


£ 2 ■ x2  y2  «*.  s2  ; 


& 


tan’ 


-1 


Vx2  +•  *2 


X = tan-1  */x. 


and  letting. 


r = r. 


- jliv*  +•  e > 


represent  the  vertical  notion  of  the  sphere,  oondition  (1.6)  becomes. 


w M 

TT  s -AO-;,  *.  cos  Q (l.fl») 

' ' f 8 ^ i 

This  auxiliary  problem  is  not  yet  oonpletely  formulated.  Since  we 
are  oonoerned  with  an  infinite  region  it  is  neoessary  io  specify  conditions  ax 
infinity.  One  oan  show  by  a simple  argument  involving  energy  conservation  that 
the  ways  motion  due  to  a bounded  obstaole  should  decay  like  (xz+  as 

(xZ4-  sz)  t*£ds  to  infinity  along  the  free  surfaee.  Moreover  the  work  of  frits 
John(2;  shews  that  this  oondition  leads  to  a oorrsotly  formulated  problem. 
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that  is,  there  is  one  and  only  one  solution.  We  thus  presoribe  the  further  con- 
dition, 

W = 0 ( 0”1/2)  as  x2+  i2— > oo  with  bounded  y (1°7) 

Eie  problem  of  the  fcraed  osoillation  of  a sphere  inToires  syswetry 
about  th*  y-axis.  It  is  oonveniert  in  numerical  work  to  make  use  of  this 
symmetry  in  the  introduction  of  the  stream  funotion  V(x,y.r).  If  r^  = x2  ■+•  a2, 

V it  related  tc  I by 

- n (loS) 

In  terse  of  V,  the  oonditicn  (1.6')  becomes, 

V - -L  a cr  7^  * cos  E©  (1.6» ») 

Onoc  the  forced  osoillation  problem  has  been  solved,  one  may  proceed 
to  the  case  of  the  freely  floating  sphere.  ?he  relooity  potential  of  the 
inoident  plane  ware  is 

W<41  = A^"K^u"Ky  (lo9 ) 


where  i is  & constant.  Ws  assume  the  motion  of  the  sphere  tc  be  of  the  form, 

y * £ j*  . where  Y,  and  pc  are  to  be  determined.  We  denote  by 

the  solution  of  the  corresponding  foroed  osoillation  problem  and  suppose 
the  total  velocity  potential  to  be  given  by. 


W - w 


on 


+ vs 


if) 


+ vv 


(')') 


is  ohosen  so  as  to  satisfy. 


J wa> 

* f 


o r» 


p - ^ 


[1,10) 


Tram,  the  resulting  solution  Y,  a nd  « nay  be  determined  from  the  equations  of 
motion. 


The  results  of  our  analysis  will  be  a determination  of  a wave  funotion, 

W,  satisfying  the  oorreot  boundary  oonditlons  and  behaving  properly  at  infinity. 
The  question  of  uniqueness,  however,  remains  partially  open,  sinoe  the  work  of 
Frits  John  assures  a unique  solution  to  the  problem  of  the  freely  floating  sphere 
only  under  the  oondition  <■  3/2 


2.  Formal  8olutlony 


▲ diseussion  of  the 
is  given  in  Appendix  I.  It 
ing  equations  (1.4),  (1.5), 
oharaoteristio  functions 


functions  involved  in  the  solution  to  our  problem 
is  shown  that  in  order  to  find  a potential  satisfy- 
(1.6")  and  (1.7)  one  is  led  to  a superposition  of 


we*  - 


9 * 

V 


■<  Pe*_,  f t * s 9) 

z.  m pi 


■ft  = 


2, 


(2.1) 


4. 


where  the  I”e  are  Legendre  polynomials,  and  the  potential  S W*  , due  to  a 
point  souree  of  strength,  5,  at  the  origin. 

W*«  4 K.(>0  (2.2) 


where 


H6  ( k * ) s Ja  ( k r\  'i  + x j,  «'k  a ) ..  c 

The  corresponding  stream  functions  are. 


***>*>«  V 


v.»  =t  ^2 a (coi  dj  — cos  3 (cos  & 1 k ^t*^coss)-  cose  Q&l* i 


(cos  e) 


2A 


(2.3) 


v*  = - •*-  2 tcxTTA  i/'Vica)  •+  !t*L  f 

P ’ -ir  J 


k .»  i n r^j  »t  cat  r-jj, 

T *"  +-  K 


a k^t'^r  nt  a*  t- 


(2.4) 


-a  »•' 


with  k,  ca-A  > „ For  the  stream  f'motion,  Y,of  the  forced 

osoill&tion  problem  wo  take. 


% ’ 


21  *£*  vz<w  <p,®) 


2 m. 


(2  o5) 


The  constant,  8,  oan  be  related  to  the  wave  amplitude  at  in- 

finity. From  equation  (1.1)  together  with  the  asymptotio  expression  for  the 
Eankel  functions,  we  find  for  tha  wares  produced  by  the  source,  S w#*  , 

-n. « i-fr  Vj-^va. 


or  writing  b - Va  V.m  , 
S - 


t±/Z 

p <r  V ? -rr 


£r  v Eir  * 

Substituting  the  series  (2.5)  into  equation  (1.6,,}  gives, 

+% 

F CO 5 2(5  = V0*fA,®)  ♦ Z ve~'*'9) 

I 

where  F is  independent  of  © . F is  determined  by  setting  <?  *■  , 

Fa-  V.Va*41  + £ _£2_  P£^(„, 

and  may  be  eliminated  between  (2.7)  and  (208)  to  give, 

mm 

V.*f.l®)e.vf<al»4)c.,  25  * 21 

I 

W' mi > » - f » s » P? ft. . .. .)■ - *<« »)  P • ■ * *„£ f »■ »>  - t-tj <. , « > - i <•  <, s 4 SJ } 


(2<>S) 


(2.7) 


(2.8) 


(2.9) 


The  procedure  now  is  to  determine  any  desired  number  of  the  by 

a numerical  process.  Such  a process  involving  solution  of  linear  equations  by 
an  interation  process  is  oontained  in  the  oonstruotion  used  in  the  next  seotion  to 


5 


verify  tba  ooavergenoe  of  the  solution*  However,  & more  oonvenieat  process 
is  that  of  least  squares  sinoe  it  leads  to  symmetrio  matrioes  with  a cor- 
responding decrease  in  the  numerioal  work*  It  might  be  remarked  that  the 
rather  cumbersome  nature  of  the  -Fg^'a  is  a oonsequenoe  of  the  faot  that  we  do 
not  know,  a priori,  the  value  of  S and  it  must  be  eliminated  in  the  computa- 
tions* 


's  are  determined,  we  oan  find  P from  Equation  (2*7), 
and  oomparison  with” equation  (1.6”)  yields. 


ODoe  the  A 2/y*  •« 


e r 


For  the  motion  of  the  sphere,  we  have 


T - 


r <x~ 


(2*10) 


(2.11) 


Making  use  of  Equation  « *2)  tu  deieraiuo  the  pressure  we  oan  also 
ooniput*  the  vertical  foroe  on  the  sphere  in  the  fora. 


R,/*')  « 2a> 


f — A.&T  C r 

I fcfa.*,*)  * “i  | T*-~  ftlJ  Wfa,8)  JikS  co$«  c/® 

A A»/v  so 

* “ 3 <r  ? a'a  J .( w/'ra.e*  + 21  cos  8 s i * a 

= - *r<|  faw;a1«)c»ssiift«  + 51  — 2—~r 77i^w} 

'“4.  V con*-  / 


(2.12) 


(2*15) 


a.  - 2 a f g.  «*  S a 

Comparing  equations  (2*11)  and  (2.12)  we  may  rewrite  the  foroe  ae, 

F„<*W  As.**/  1*2*1  ^ **<**/  fMDj  rfjr 

v'q  l ,F|*  j siA*  1 iFl*  7 dx 

with  bars  denoting  oomplex  eon jugate s*  This  latter  fora  is  <he  familiar 
decomposition  of  the  foroe  into  a term  involving  a virtual  mass  times  ac- 
celeration, and  a damping  faotor  times  velocity* 

Having  oonsidered  the  foroed  osoillation  problem,  we  may  turn  to  the 
motion  of  a freely  floating  sphere.  We  note  that  in  the  above  solution 

m«  m #4  «— .*%  A V«r*  ~ V- V • WT^  r -I  IV  . J JLt-  - . - a 4.1-  - -I  4 - 

**-*•-**  w’j»  m mmwaawu^o  w a same  ymxuww  wo  ^4  xuou  uaa«9  ^uao^  ui  UlO  Nava  MVWIViA 

at  infinity).  In  order  to  obtain  the  potential  for  an  arbitrary  phase,  <*  , 
which  is  needed  in  order  to  prooeed  as  outlined  in  Seotion  1,  we  replaoe  t by 
t + in  the  previous  result*  *he  corresponding  W will  be  denoted  by 

and, 

WS  * w' 


where  W is  the  solution  obtained  above* 
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Our  task,  then,  ia  to  find  a funotion,  W**'  , satisfying  equations 
(1*4)  and (1.6)  as  well  as  equation  (1.6)  with  at^  =.  - »®tion 

is  no  longer  sysmetrio  with  respeot  to  JL  and  requires  nore  general  functions. 


( P?~<cct6)  A 
( pia»*  £*-*.»*  (*24  )*  **  '>*<*  ■*  j 


= ft-; 


-i  / f P*T  Vcssfi) 


vg-'.'SSSaJZ: 

«a  l f>2++g 


2A-2'K-‘2.  ~P*A**  ^ 


^ 3 «,2,*  - * *,  "*  = i.Z,  • • A 


the  P^(cos^)  being  associated  Legendre  functions. 

With  y t (o}0  , x*  (p  sin  3 cos  jt  , we  have  from  equation  (1.9) 

>wca 

L-i*  ~K^X  CCS*[[cos5  (ts^a«inee(>a)  + smScojX»  m(Ko.jir\afoiX)] 

4-  x fcosfi  i m(K4  sitidfosX)  - sir\6  cos^c»{  («  a s ini'  cos  X*)J  j 


=•  - k A a 


- ko.  c o s 9 


{ ftfKd,  a, pc)  + -a  .T  ( k o.,  S,  x)  } 


How  £ and  I are  even  functions  of  /.  with  respeot  to  ?C  - 0,  hence  each 
will  have  an  expansion  in  { with  j&  * * x replaced  by  oos  . 

Moreover,  R and  7 are  even  and  odd  respectively  with  respect  to  X * "Vk 
henoe  we  take. 


w“".KA  I I - -2 


■ * k a i f 

***  I sns  l I ^ ft. 


* (*£'«  j 


CSS  C/v«X 


2<w-2'm+2  p 


£iifiS4Llcos^-,W 

pt*+l  J 


with 


(2.16; 


ISJ^jJT.ttolSl}  c»»  i-w*  (2.17) 

A-‘—X  £<>..*,  C/»»«-gSgS^  C-V...)}  ...  (2.15) 

The  dependence  on  ?t.  in  equations  (2.1?)  and  (2.18)  aay  be  eliminated 
by  asking  use  of  the  orthogonality  of  the  oos  m\  pc  to  write, 

AA  «♦**■  E'n  - z*v\ 


4,.,-Hfccos®  r* 

- it  , J0  *•?*?*>  c°*  C * . /4£/^  < K * > 


(2.19) 
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Z *\Z  #£^(Vcosd> 

■'*.  * '/•■A. 


2/v*  - 2/vt*  ♦■'2 


- -K4C6S  S 
5 — X 


r*'x 

J I(K4  jCjX)  C6J 


(,2/v*v-i^«7C  = ©V 


(2.20) 


for  s " 0,  1,  2,  ....  ill  (2.19),  1,  2,  3,  ....  in  (2.£o)  Finally  we  may 

make  u«s  of  the  identities. 


-Tjn 

I C»$  Co3^\  C<)$  2 '**'  X - ~ C - i)  >VV"'  Cts-« 

^ 2. 

//? 


| 5in  c « S xi  CoS  te-m  - 1 1%^x  = ^ i’- ( ( ») 
to  carry  cut  the  integrations  for  the  functions  kM(  tea}  & ") 

A ^ s f«  i)  jl  ^ coif  f s a f (n  6 I + S • h (9  /Xu  ( ^y<cn,  tir\  Jl  * 

At^_*  (' '^’K*casfl[£  co?£  ^/wsino^  J-Sir»5  [ sme>- 


■ J**-(**  *•*•0}  (2.21) 
(2.22) 


The  determination  of  Vv***  is  thus  reduced  to  a numerical  process  which  is 
analogous  to  that  used  in  finding  v/*5  • Once  is  computed  the  vertical  foroe 
on  the  sphere  may  be  found  from. 


W 2H 

The  dynamio  equation  of  motion. 


(2.23) 


Tr 


3 + t ra  so .*  ? 


(2.24) 


will  then  yield  two  equations,  by  equating  real  and  Imaginary  parts,  for  the 
determination  of  "fc  and  u • 


3.  Convergence  of  the  fbrwal  Solutions t 

The  aim  of  this  section  is  to  obtain  some  results  concerning  the  converg- 
ence of  the  inflate  series ’used  to  construct  solutions  in  8eotion  2.  The  analy- 
.1.  "ill  ho  iuaIo  ul  Lima  Laly  Lu  i cel  ou  wall  know  a resulto  in  the  theory  or  ex- 
pansions in  series  of  Legendre  functions,  however  certain  manipulations  must  be 

performed  before  one  is  in  a position  to  invoke  these  results.  We  will  deal  only 
with  the  foroed  osoillation  problem.  Although  similar  procedures  should  carry 
through  for  the  freely  floating  sphere. 


Before  proceeding  we  wish  to  make  plausible  the  type  of  expansion  that 
was  adopted-  In  Seotion  2.  For  this  purpose  we  consider  the  two  simplified  problems 
which  arJ.se  when  cse  lots  o~  become  infinitely  large  or  small.  The  free  surfaoe 
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condition  (1.6)  may  then  be  replaced  by  the  simpler  conditions  W s 0,  and  Wy  * 0C 
respectively.  In  addition  the  conditions  (1*7)  nay  be  deleted  and  it  suffioes  to 
demand  that  ¥ be  regular  (i.e.,  vanishing  like  l/p  ) at  infinity.  Consideration 
of  formula  (1.6)  of  Appendix  I shows  that  w£  — * 0 as  — * «»o  , and  W0*-*  ^ 
as  <r  — * 0.  Moreover,  we  have. 


K 


W, 


P,  ( e o 5 C ) 
(° 


as  <r  —*  oo 


(3.1) 


w, 


tM 


(£  • t &) 
p tit  4 i 


2 % 


V —4  O 


(5.2) 


Vow  it  is  possible  to  solve  the  two  limiting  problems  by  standard  methods 
of  separation  of  variables.  For  the  case  f ■ 0 on  the  free  surface,  the  solution 
involves  a constant  times  the  dipole  potential,  , which  will  fit  the 

limiting  farm  of  our  expansion.  Similarly  for  the  oase  Wy  s 0,  we  oan  fit  any 
given  funotion  over  the  immersed  surface  of  the  sphere  witn  an  infinite  series  of 
even  indexed  Legendre  polynomials.  Hera  again  the  leading  term  is  provided  by  W0* 
and  the  further  terms  by  W 2 ' s . 


These  remarks  indio&te  also  how  the  speed  of  convergence  should  depend 
on  v*  • ' For  the  cars  of  large  <r  , the  leading  term  above  should  be  nearly 
adequate.  If,  on  the  ether  hand  the  calculations  for  the  oase  <r •*  0 are  carried 
out,  it  is  found  that  the  coefficients  of  the  PJ/v^ :otSj  are  of  order  n"5/2  eo 
that  the  convergence  is  slow.  (Xxr  method  will  yield  a proof  of  oonvergenoe  only 
under  the  assumption  ox  small  o~  which,  however,  from  the  above  discussion  will 
be  the  oase  likely  to  give  difficulty. 


A precise  statement  of  the  main  goal  of  this  seotion  is  the  followingt 
Theorem;  If  £ A «osd  where  k is  a oonstant,  there  exist  constants  S,  A£,  A4,  „» 
such  that  the  series, 

21  * 5 { W*  +•  £ a.  2”'* 2 y 


has  the  property 

. d w**  ** 

l Tflp-"  Z [(Z« +*$„<'■<>**)  - ka 


til 

<>f  I 


(»  t c. 


(3,3) 


the  oonvergenoe  of  the  infinite  series  being  uniform  and  absolute  in  c « **  -n/,,  * 

Suppose  first  that  equation  (3.3)  holds.  We  obtain,  thenr  by  setting 

© « 0, 


■J"  — * ^ " 21  0 PffA  ^ 1 * ~ - 1 ( 1 ' U 

Sinoe  Fn(l)  • 1 for  all  n we  obtain  by  subtracting  equations  (3. a)  ana  (0.4). 


(3.4) 


• -7-* — - +- 


0» 

- Z {(Cm-ulf  P2<^c«5#)-cos^^  ■**(^5,*  * ® 1 ~c  0 5 •»}  (3.6) 


We  should  remark  at  this  point  that  the  series  (2.7)  is  obtained  from  (6*5)  by 
integration  so  that  the  uniform  oonvergenoe  of  the  latter  implies  oonvergenoe  of 
the  former. 
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Now  the  left  hand  member  of  equation  (3*5)  ie  a function,  continuous 
in  0 $ £ < -nv£  t and  vanishing  on  9 « 0.  As  such  it  oanbe  extended  across 

® /2  as  an  even  function  of  S , to  define  a new  function,  continuous  with 

a piece-wise  continuous  first  derivative  in  o £ 5-  a.  vr  . it  will  have,  then, 
an  wXpMMSivM  of  tbs  fora. 

*L  ae-wv  + ,)  “ Ofwn-^) 

^y\  - 4 

and, 

09 

a#  “ - 2T  (3.6) 

If  the  series  on  the  right  converges  uniformly,  it  defines  a funotion  with  the  same 
properties,  and  hence  will  possess  a similar  expansion, 

ZT  '*'***  1 C 0 S 9 ^ 

O^SO 

where  the  's  satisfy  (3.6).  If  equality  holds  we  have  necessarily 

for  a * 1,  2,  ...  and  therefore  also  a0  = b-. 


The  proof  of  our  theorem  consists  in  giving  a method  of  construction 
of  a sequence  { A^-}  for  which  the  series  does  converge,  and  -tr^s  «Iwfor  m ::  1,  2.. 


We  make  use  of  the  formula, 

/>”/* 

J Pa*»€*S®)  cos  6)  do 


if  <w  * *\ 1 j 


- o if  n/n'  and  n-n»  is  even 


/*  01*1, 


if  n - n*  is  odd. 
n even,  n*  odd. 


Then  if  equation  (3.6)  holds. 


^2/w\  ~ — 

4m  t-  l 


2»  : *■ • 


’{» 


+ Ka  < — i V 


e^em-ore 
a -m  ! r e<*>  - o • 


1 r* 

♦ 2 Km 


A C2/*\  f i - ka). 

Z/V\ 


i 


fen»\  + + fa  K'm  1)*  - •)!  8 


* h 2 -• 


(3.7) 


and  our  problem  is  to  establish  the  existence  of  a solution  of  these  equations. 

In  order  to  prooeed  we  Interpret  equations  (3.7)  as  a single  vector  equa- 
tion. Denote  by  x ° fljsd  X the  infinite  dimensional  vaotors  having  and 

as  components,  i.e.,  x«  = -)  , *,£  c a.^  j xrC*„Xj--. 

A.  L2l±J  Further,  by  P f?)  denote  the  matrix, 

S'*  'lAv  ^ I 

p _ S'V*'  1 + “KS.)  (-I)*"'*-' 

2 —>)  Czrr*\  + z ) f/Yv\  (z  m t 

(,4/n-t-l  > K cc  A/m  * tt/*- O 

*7  w — i ■ — — — — -—  ■■■■■■■■ 

^ Jm-f  Z'Uft)  (Z«rr*  - 5,yi  *.  (jf  Z/w  +Z*V>  K'Sw  I2,  ( r*'~  ' ^ • 2 


1C 


and  than  the  equations  (3.7)  beoome. 


X - x°  +•  Frx) 


(3.7«) 


Consider  now  the  set  of  all  infinite  dimensional  vectors,  x , 
satisfying  the  condition. 

Oft 

U X fl  — ^ + l j I Xm  1 K.  / 2 0 J 

■>>■=■  l 

This  set  forms  a linear,  normed  vector  space  I,  whioh,  with  respeot  to  the  norm 
defined  in  (3.8),  is  complete,  i.e.,  is  a Banach  spaoe*  % this  we  mean  that  if 

f x*  V is  a sequenoe  of  vectors  of  X,  satisfying  nx"-x'w,'Ho  as 

there  exists  a vector  X , also  in  X,  suoh  that  * -xMti  = o.  Prom 

equation  (3.5)  the  vector  x°  defined  by  the  's  is  contained  in  X,  and  we 

prooeed  to  show  that  equation  (3.71)  has  a solution 

Lema.  t If  B is  a Banaoh  space,  and  F is  an  operator  on  B satisfying  a Lips  chit* 
condition. 


II F ( - p l li  4 M « > - ^ « ocr-i<i  (3.9) 

then  the  equation  (3.7* ),  for  ”*•  e B , has  a unique  solution  Sea  . 

Rroofi  refine  the  sequenoe  {x^y  as  follows. 


X ° 

, X 1 = x • F(X»)t  A z 

: X'+  Ff*'l,  ••  • X~+'s  X»  4-  FCxm  ) , • 

y ' 

- X • - 

F ( x*  V 

x* 

- X.'  •= 

FCZ')  - PtJM 

/,  BX*-  X'll  < M frFfx*)ll 

X 3 — X1  s 

Ff^-Fr;'| 

e 

II  X * - X*l|  £ M II  -X  M|  £ M*  !• 

k * 

x ~ 

ft 

♦ 

- X~-' 

• 

ft 

- F ( ) - F(  *•*-*) 

• 

e 

• 4 
n x y 'v'”  ‘ il  < m ,n” 1 n rrxM  U 

'For  <y\  > iw\ 

X ~ ■=■  CFPxm"')-  - Pfx  *’*)] 

+ * . . + C Ffx  - Ff 

/.  I|  x~_  x~i  <,  M II  X «*-•-  x ~*Mi  + M ft**"*’-  x «■  * ' * M H*  2 **-'  1} 
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S ‘ 4 M""]  « ^-+M4  * • - -* 

4 M'"  F^*)  -~  ^ 

• • HU  **  X <w'  II  * ) 0 2 S <va  , **“  V eo 


Hence  there  exists  z suoh  that  ti  x'm— x n — * o sis  s i — * * Moreover 

from  the  Lipschit?  condition  we  immediately  infer  F fit  ) +•  x * = x 
well  as  the  uniqueness  of  the  solution. 


In  order  to  enter  the  form  of  the  lemma  we  must  verify  that  the 
matrix  satisfier  the  Lipsohitz  condition.  We  have. 


FC *8 )li 


* z z 


“2  AV* 


£ ^ ci~~->ne+» 


VS 


( t ♦ K o-)  0 I X/k  “ l 


as 


+ KG. 


ST  y g/»A  i g m ■ (H/n  + O / --sj„ 

^ 2*m  + *‘'**(2/**  ♦ >)  f£<W> ♦ it  i )Z  (e»  i ) C'V  “ I ^ * 

ee s i 


(3.10) 


We  make  use  of  the  numerical  series. 


2-*  J 


l Z2~  (/*  1 }2  2<v\*.( 

and  estimate  as  follows \ 

r-  2 '*•'  1 


i v 

2~  1 


. . . «t 

• $ * fc.rw\  ▼ a i »•»"  # i 


2 («*-')  \ 


‘ £ £ 


E <<***  ! >*■  tz>wi+  0 


- -i  ( 2: 

4 1 l 1 


y a v*  1 <t  2.  T a ** 

“V,  ' 2-  , 


f 2/r**  i ) 


- { HL-o. 


From  (3.10)  we  obtain,  then. 


Hr^r.  > • r*  / — . s » i / rv 

**•  v * ' • ‘ a ‘ “ ^ 


t 4 


4/  a I il  w ,«.  . M 

j ••  *'  a •• 


It  follows  that  the  Lipsohitz  oondition  (5.9)  is  satisfied  fors 


K a.  < 


f r 


_ JL  * 


f - i * 1 

or  approximately,  K ^ < 2 


(3.11) 
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We  oonolude,  therefore,  that  the  equations  (3.7)  possess  & solution, 

Ag,  A4, . • • such  that 

5?  * (AttA 4,  • ■ % ) « X , v*.  s O (*,.“*-*)  , € >e 

It  follows  that  the  series  in  the  right  hand  member  of  equation  (3*4)  converges 
and  we  use  that  equation  to  define  the  constant  S.  Moreover  the  « ’s 

haw  been  determined  so  that  equation  (3.5)  holds  and  by  multiplying  (C-.i)  uy 
eos  W and  adding  to  (3.5)  we  dsduoe  that  equation  (3.3)  holds.  Ihia  completes  the 
proof  of  the  theorem,  for  suf fioiently  small  < *.  . 

r 

A solution  to  the  problem  of  a freely  floating  semi-immersed  sphere  is 
obtained  in  a form  which  admits  of  numerical  oaloulationB.  The  solution  con- 
sists of  fitting  an  infinite  series  of  non-orthogcnal  functions  to  a given 
boundary  condition  on  a portion  of  the  sphere.  Some  estimates  are  made  regarding 
convergence  and  some  affirmative  results  in  this  direction  are  obtained. 

Once  the  numerical  calculations  of  the  ooeffioients  ere  carried  cut 
certain  quantities  of  physioal  interest  may  be  obtained,  e.g. ,the  wave  motion 
at  large  distanoes  from  a sphere  moving  with  foroed  sinusoidal  oscillations,  and 
the  virtual  mass  and  damping  factors  for  a freely  floating  sphere. 

As  in  the  previous  report  on  the  oylinder,  oertain  approximate  results 
could  ts  obtained  tot  the  sphere.  One  could,  for  example,  solve  the  simpler 
problems  in  whioh  t'ne  free  surface  condition,  is  replaoed  by  W s Wy  s 0 
corresponding  to  large  and  small  frequencies.  Also  a combination  technique  oould  be 
used  in  which  one  negleots  the  diffraction  effeot  for  small  frequencies  using  only 
the  potential  of  the  foroed  osoillation  problem,  and  the  simplified  condition  W * 0 
for  large  frequencies. 

The  simplified  problems  can  be  solved  by  Beams  of  Green’s  functions  vhieh 
oan  be  expressed  as  source  systems.  We  have  cot  included  these  results  for 
two  reasons.  First,  the  eolution  obtained  in  this  report  seems  to  be  complete  in 
itself  without  the  neoessity  of  appealing  to  approximations.  Seoond,  these 
problems  lend  themselves  to  variational  techniques  which  enable  one  to  determine 
the  changes  in  the  Green’s  functions  with  small  changes  in  the  shape  of  the  obstacle. 
It  is  the  author’s  hope  to  disouss  this  approaoh  to  the  problem  elsewhere. 
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APPENDIX  I . 
HAVE  FUNCTIONS . 

1.  Characteristic  Functions i 


We  are  oonoerned  with  the  determination  of  oertain  functions.  vrhioh  we  oall 
wave  functions,  which  will  permit  the  oonstruotion  of  a velooity  potential  satis- 
fying equations  (1.4),  (1.5),  (1.6)  and  (1.7).  As  a first  step  we  indtoduoe  char- 
acteristic functions  whioh  satisfy  the  first  two  of  these  equations.  These  arise 
in  a natural  way  if  one  Attempts  a solution  by  separation  of  variables.  Consider 
the  functions. 


wL"w)  * { 7* ““l1  - V*“  ^ 

elA*'  J • * - 


^TVasxl  w f - k ....  puC('cosS>)7 

r I ctA*t  j"'4  **ai*V**  ^-<>2., * 

where  P£Vcoy<4)  are  associated  Legendre  functions.  These  functions  are  harmonic 
in  f > o'  , and  we  proceed  to  show  that  they  satisfy  the  free  surface  condition 
(1.6) 


C 

• t 


**  *•  »;  !j  ’ ’ • k - 1 

(r . n 


* 


'V*T 


is  xei-o  if  n m is  odd.  Similarly  from  the  reoursion  formula 


so  that  P2T  c 0 ) 

dlff  I ■*-(/* ,+  ,)  P~to) 

**  U.  _ 

dL  P — . 

we  infer  that  — — — is  0 if  a - n is  even.  Using  these  results 

A x 


±±ij£t 

* <)  a 


9:±VX 


‘ vff/  • f 

'•  *1% 


2-A,  “S'ws 

2 A - 


pIT  t*) 

— -*■  =6 


lAti 


4 - 


l'  * * 

f-  K . K } 

■ < 


It  is  apparent  that  a series  of  the  ^1**  will  never  be  adequate  for  our 
problem  sinoe  such  a series  would  vanish  at  least  liks  Q "2  for  large  x»  , 
oontradloting  the  oondition  (1.7).  To  eliminate  this  difficulty,  we  mtroduoe  a 
"source*  solution  whioh  will  satisfy  equations  *1.4)  and  (1.5)  and  also  behave 
properly  at  infinity. 

2.  Source  Solution t 


To  obtain  the  desired  solution  we  solve  equation  (1.4)  when  a point  souroe 
is  admitted  at  ( x*,  y',  iJ),  that  is,  we  soek  W*  (x,  y,  z,  x-,  y1,  *») 
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satisfying. 


< - J <*-*')  U*j  -^*)  SCi  - *')  (Io2) 

whore  the  right  hand  side  represents  a produet  of  delta  functions.  Ve  assume 
initially  that  the  souree  point  is  (o,  y',  o)  so  that  W*  will  he  independent  or 
6 • Equation  (l.2)  baeoaes. 


a*w« 


S f o'i 


Applying  the  Henkel  transform, 


- j"  d.  Jp  (n  J,)  W^Ca, ^ ) oIa, 
yields,  ° 

Tif*  " k2jAr  * ~ - - *'-*  *'> 


(1.3) 


Equation  (1*3)  is  an  ordinary  differential  equation  for  w(p,y)  wnion  may  be  solved 
by  application  of  Fourier  transforms*  The  resulting  solution  is 


The  boundedness  of  f as  y — * « and  the  boundary  condition  (1*5),  yield. 


jwj-  C 


= T f 


► “ * 


*1 


(I.*) 


The  transform  w(p,y)  is  regular  in  - <%  except  for  a 

simple  pole  at  p * K.  The  ohoioe  of  a contour  for  the  inversion  integral  is 
dictated  by  the  behavior  we  want  at  infinity*  We  choose  a contour,  o,  which  con- 
sists of  the  positive  real  axis  exoept  for  a small  semi-oircle  in  the  lower  half- 
plane, about  the  point  p " K.  Tnen, 


* J[  ZCf,*)  | jL£JL  jb- I*  + j d l.  (l*5) 

Using  the  identity. 


J T^JWU^*  a 


4 

VZfiTjr*- 


and  replacing  n.  1 by  obtain  finally, 

v ^ o » - IX 

C 

R*  =•  (*-><•)*  +•  n1  + (*  - t1)1 


where 


(1-6) 


The  function  W*  is  thns  seen  to  be  identical  with  the  Green's  function  ob- 
tained by  Fbrits  John(^).  It  ia  interesting  to  note,  in  passing,  that  it  nay  be 
found  by  a method  used  by  Xennard(4)  in  the  case  of  two-dimensional  motion*  In 
this  procedure  one  first  considers  -the  simplified  boundary  condition,  W • 0 on 
free  surface,  y * 0,  and  treats  the  resulting  problem  as  an  initial  elevation. 
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It  08ja  be  seen  by  direct  differentiation  that  W*  satisfied  equation  (1.4) 
and  (1.5).  In  order  to  demonstrate  that  it  also  satisfied  condition  (1.7)  we  mak6 
a transformation  which  at  the  same  time  yields  a form  more  amenable  to  numerical 
computations . From  equation  (1*6)  we  can  write. 


= ±.  * J/KA)  R*  PaV  f L Z2  4*'*^  Uj  I f>  a ) A |e> 


where  h!)'  is  the  Hankel  function  of  the  first  kind.  Again, 


W„  -L  - e^KJt 


c‘  P 

where  c is  the  image  of  C in  the  real  axis.  Sires  O ( ^ 


VT 


for  large  >%|  , we  can  more  the  oontour  of  integration  to  the  positive  im- 

aginary axis  and  obtain, 


W*  r — +iir  « * 
i\ 


w*-*'),.1*!  \ 2 flrrx-  k*  - 1 K w). 

^ ^ nt  oca)  *T|  1 '-K.ft.Wr 


where  Koi:a)*  is  a quantity  which  is  real  for  real  t , Finally 

wo  use  the  identity. 


J K 0 ( a.  Jt~)  cos  -6-  a el*  — .X 


VaJL  + j,'1 


/ w n \ 

v±.  I ; 


to  write. 


W*=  +■  ^7  -»?ttka  H^Cwa)  4k  r icc»s^^+qO  + 7 

* o Tv+K'- 

where  A,{  : ;A.X')1  -h<y-y')z  t*  <*  - * ')*,  Now  = o<V* 

and  henoe  we  see  from  (1.8)  that. 


(i-8) 

) for  large  r. 


- 2 je.  * + 1 * -A.  K n 4.  o C -4 ) as  a — * <w  y bounded. 

a r ' 

Some  further  steps  toward  numerioal  computations  are  given  in  Appendix  II 
where  it  is  also  shown  that  the  apparent  singularity  in  equation  (1.8)  vanishes  by 
subtraction. 


5.  Axially  Symmetric  Motion: 

Suppose  the  motion  is  independent  of  X.  It  is  then  permissible  to 
introduce  the  stream  function,  defined  by 

dy  ^ n.  ^ 

3a  3a 

The  stream  function  v **  corresponding  to  aa/ * , 


(1.9) 


with  x‘*  o is  given  by. 
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^ - ^ri'1  + P tUL  A-^+v>A  oij> 

r " * 


(I. 10) 


as  may  be  verified.  on  noting  that. 


d\^  T,  ' A.  jj  T,^a)  j Te  ( |>a)  = - £ or,  r p 

me  integral  term  in  equation  (I. 10)  may  be  manipulated  as  before  to  yield. 


tZlL  - JL  f « i r» 

+ "Tr  ® 


aj*)  A.  K,  fTA)  «<  T 


'T  "*■.+■  K 1 


.A  K,fTA)<<T 


■where  - -j£  M*,'1  <<*■*:  '».  . The  first  integral  may  be  evaluated  by 

noting  that,  A.  k,  tta)  is  taI  , integrating  by  parts,  and  using  the  identity, 

(1.7),  quotad  before^we  obtain. 


•v*-=  - 


< 'H  -V) 


~i  — - - 


+•  ZVK  X 


40 

4k  f k «•  a 'r''jj  + Y*  - T c o a r<'-a  + V* 

T 4 T*-  *-  XT  * ^ ! 


/)  *tnr 


(1.11) 


Since  the  motion  is  to  be  independent  of  JL  the  only  charaoteristio  functions 
to  be  considered  are, 

f y,L<co*S)  K r£l.,uosfl'  ■) 

- TT  ““-ip* S **<*«.••• 

The  associated  stream  functions  are. 


y-  - pu-ftl,5fi1  ’ ces9  P£*Ccos©\  k P^dosS)  - CotO  Pn-fc  o$6\  ^I#X2) 

frs.  a*., 

To  see  this  we  note  that  in  spherical  oo-ordinates  the  equation  (1.9)  reduces  to. 


i 3v  _ „ ^ <?W 

o r-»  — ~ P t ^ — — 

8 ^ V Ap 


and  using  the  identities. 


17 


(m- mI  R*  ( 


d X fit* 


i^.- 

/X 


it  may  be  verified  that  s/Sx  and  Wj.^  satisfy  these  equations. 

The  introduction  of  the  stream  funotion  V proves  a convenience  in  solving 
the  forced  oscillation  problem  of  Seotion  2,  since  it  eliminates  the  neoessity 
of  numerically  computing  derivatives  of  TY*.  The  results  of  the  last  three 
sections  lead  one  to  choose  as  the  stream  function  for  this  problem. 


V = s V*  -K  Z_ 

I 

where  V0*  is  the  stream  function  for  a source  at  the  origin.  The  origin  is  again 
chosen  as  a numerical  convenience  since  it  3l:-.ghtly  simplifies  the  computations. 


APPENDIX  II 

AIDS  TO  NUMERICAL  COMPUTATIONS 


Numerical  work  based  on  this  report  requires  computations  with  the  two 
functions. 


w»%,a  - f - te  ( — (n.i) 

v,v.»  - c“-23  (n.*> 


In  this  form  the  calculations  are  lengthy  beoause  of  the  occurrence  of  the  trig- 
onometric functions  in  the  integrand  and  the  infinite  range  of  integration.  By 
suitable  manipulations  it  is  possible  to  reduce  these  expressions  so  that  they 
involve  only  finite  intergals  with  monotonic  integrands. 

For  wo* f p - we  need  the  two  integrals. 


R •=. 

wVana 

••  **■>**  v 

R = 


1 f tiiJL*!  K0cc-al  d T 

<4  v--*-  K v 

jii  . ^ i 2 i f.  5 

w 

r c#<  x-w 

' I — s— , 

om 

I A*.  S ( r\ 

J - — K (Ash'S 

>4  * a *•  * v 


9* 

s - f ^JLL"JL3  K#  tk  "t* 

'd  Tl  f K» 


«.  a.  c o c w n - 
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Now 

Ke(/j*0  « f — ol* 
A • vT~ 

end  we  can  write, 

R A [*  f c # 4 4’(£j*' 

* [ VJT^  l ^ — 


+•  t o s ^ 


l »t-pk.\ 

--Vot  ~ 


p-«b  <•  •< 
> K 


/4.  <■  f 1 


oo  *o 

s 1 f oU  f Sin 

"J.vPrrl  — 


A1  + l l 


{ I >\  « ( k ifc  ^ 

s>n  <**  fft*:  -*  1 j.. 


(o*~  <•  •< 

P*  >-< 


Hromth6  identities. 


W 

L 

r 


* ^in/AX  j ^ _ jr  ^ V 

xi+-<xi  2- 


C_LL^_X  ^ 


for1  > i.  , o*  ^ o 


’fa  i*  __/iA  > o , ft  >,  O 


it  f ollow8  that. 


90 

/v»  p - * *«* 

tSU+SIl  •*■ rf* 

'0  VtTTT 


t Z K i 


^ a/^”»  T 7 


•ot 


Moreover, 


f ±1  *— b ? f $,<**}  - V**'] 

jk  A 

* f I 

where  5^(4)  is  the  Struve  function  of  order  0.  Henoe  we  obtain. 


W*l  s £ -tka 

°V:k  e 


-KfcCoS £ 


[ Setx&  $ 1 r\»)  + T.  <**  5|f'c,l 


_ .|UC0f$  f J Kft*  !"M  1 , 

2 K A ] ~ •<  * 

'<  V A.  * T-  i 


+ ffKXi 


>RCL  (Of  9 


3j  f Ks.  • > n *) 


(II. 3) 
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For  small  x. 


~ JU-^  ^ where  y"  is  Euler's  oonstant.  In 

addition  the  integral  in  (II.3)  oeases  to  oonverge  as  a — 9 o • Hence 
the  formula  ceases  to  be  of  use  at  a =0.  It  is  most  convenient  to  study 
W0*  for  small  d-  , by  means  of  the  original  equation  (ll.l).  If  we  replace 

ic0fT'O  by  its  asymptotic  estimate,  (***^)  ant*  “S^®  use  °**  ^® 

formulae, 

/ 'xTl'a'g  c 0 ~ * t & (y**-  Sc  </*•*] 


A & 


A 


*4 

JT  iiSLI 

where  £c  ( i>  \ , E.  f Y ) *r«  ®xP°a®i1tial  integral  functions  (see  Jahnke 

and  Emde^),  we  obtain. 


- 4* 


If  f K c o s + <r  s . n % . - ,<•-*  « 

TT  J ^W-TTi ° *c/T*  I T ^ 4 K Jt  * ^ PA 


2 fc  £ f ^ • 2 K Xir,  i 

a VI  Y 


( tt.a> 


The 


limiting  value  of  Wc*/  is,  therefore, 

P - 


2 -»  - KCX  O * , \ 

- 2 * 4 fcl  I'Kixl 


Similar  calculations  can  be  carried  out  on  the  integral  in  equation 
(ll«2)  and  lead  to. 


f K J IA  t-H  - 

J.  -7~~- 


T fi«  t« 


A K.  P'TA^  4 n-  * 2 c 0 S $ - A 


) A r- 


+ vKa  «;.  8 4 


• k*  s 1 ►> 0 


[ Y,  (KASinfi)  - S,  ( K a.  « l n d ^ ■*■  ■—  J 


-A. i?  4 - xi  - «x a x*  — -a a _ — a 

’“«* « wi  10  on«  oui  uvo  1 uxio  uiuu  01  ui  uoi  x • aau 


(ll»5) 


< * 1 Tr/i  N = “ K A.  L S,  f K «0  ~ i J •»  4-W  K(^  1 

V * C a , c ^ = c> 
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